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Quantum fluctuations of the electromagnetic 
vacuum are responsible for physical efltects such 
as the Casimir force and the radiative decay of 
atoms, and set fundamental limits on the sensi- 
tivity of measurements. Entanglement between 
photons can produce correlations that result in 
a reduction of these fluctuations below the vac- 
uum level allowing measurements that surpass 
the standard quantum limit in sensitivity.^"^ Here 
we demonstrate that the radiative decay rate of 
an atom that is coupled to quadrature squeezed 
electromagnetic vacuum can be reduced below its 
natural linewidth. We observe a two-fold reduc- 
tion of the transverse radiative decay rate of a 
superconducting artiflcial atom coupled to con- 
tinuum squeezed vacuum generated by a Joseph- 
son parametric amplifier, allowing the transverse 
coherence time T2 to exceed the vacuum decay 
limit of 2Ti. We demonstrate that the measured 
radiative decay dynamics can be used to tomo- 
graphically reconstruct the Wigner distribution 
of the the itinerant squeezed state. Our results 
are the first confirmation of Gardiner's canoni- 
cal prediction® of quantum optics and open the 
door to new studies of the quantum light matter 
interaction. 

The quantization of the electromagnetic field implies 
a minimum uncertainty relation for non-commuting ob- 
servables such as photon number and phase, or the in- 
phase (/) and quadrature phase (Q) amplitudes of a 
mode of the electromagnetic field. The electromagnetic 
vacuum is a minimum uncertainty state with quantum 
fluctuations distributed equally between the two quadra- 
tures. Parametric amplifiers operating in the opticaF^^ 
and microwave domain have been used to produce 
squeezed states of the electromagnetic field, wherein the 
fluctuations in one quadrature are increased and fluctua- 
tions in the canonically conjugate quadrature are reduced 
below the vacuum level, allowing for an improvement in 
measurement sensitivity.^"^ The focus of our research, 
however, is to reveal the effects of squeezed vacuum on 
the radiative properties of an atom. In the optical do- 
main, only a few experiments have explored the squeezed 
light-atom interaction, with studies in free space^^'^^ and 
in a cavity-QED architecture.^^ Our experiment takes 
place in the microwave domain and uses a polariton 
qubit — an effective two level atom formed via the strong 
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FIG. 1: Experiment setup, a The phase space of a mode of 
the electromagnetic field is described in terms of its in-pliase 
(/) and quadrature phase (Q) components. The Gaussian 
variance of the vacuum state is shown as a dashed line, and 
a squeezed state as the green region, b A lumped-element 
Josephson parametric amplifier is used to generate squeezed 
vacuum that is coupled to the input port of a 3D transmon 
qubit via a circulator with coaxial cables, c The state of a 
two level atom may be represented on the Bloch sphere with 
angles 6 and describing the latitude and longitude respec- 
tively, d The resonant strong light-matter dipole interaction 
of the transmon circuit with the 3D cavity results in two po- 
lariton states |-|-) and |— ). The bandwidth of the squeezing 
is centered about the l^r) — >■ |— ) transition frequency and is 
large compared to the natural linewidth of the transition. 



light-matter dipole interaction between a superconduct- 
ing circuit and a microwave frequency cavity. We use 
a Josephson parametric amplifier to produce broadband 
squeezed vacuum in the modes of a transmission line that 
are resonant with the atomic transition. The architecture 
of a one-dimensional radiative environment^'^^"^^ and the 
strong coupling available in circuit-QED^*^ enable us to 
engineer the radiative decay to be solely into the modes 
of the transmission line that are occupied by squeezed 
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vacuum. Thus we are able to systematically explore the 
dynamics of the atom under the modified radiative reser- 
voir. 

The itinerant electromagnetic field generated by a de- 
generate parametric amplifier may be approximately de- 
scribed in terms of the squeezing moments, M and iV, 
that are related to the frequency correlations of the out- 
put field. In the limit of large amplifier bandwidth 
(a'f(a;)a(w')) = N5{lo - u') and (a(w)a(w')) = M5{ijj + 
lJ — 2a;o), where a, are the creation and annihilation 
operators of the output field of the amplifier, and wq is 
the center frequency of the amplifier. Squeezed states oc- 
cur when M > N, with M bounded by < N{N + 1). 
The radiative decay dynamics of an atom that couples 
to a broadband squeezed reservoir centered about the 
atomic transition frequency are governed by the optical 
Gardiner-Bloch equations,^ 



(fix) = ~j{N-M + l/2){a.,), 
{ay)^~j{N + M + l/2){ay), 
{&,) = -7(2iV + l)(a,) +7. 



(1) 



Here, a^, Oy, and, are the pseudospin operators for 
a two-level atom. As shown in Figure 1, the Q quadra- 
ture of the electromagnetic vacuum is squeezed, and a 
coherent drive along this axis induces rotations of the 
atom about the y axis of the Bloch sphere. By setting 
M, = in Eq. 1 we recover the case of radiative de- 
cay into electromagnetic vacuum, where the transverse 
coherence decays half as fast as the longitudinal coher- 
ence {T2 — 2Ti — 2/7). In contrast, the radiative de- 
cay into squeezed vacuum is characterized by timescales, 
^ Ti/{N - M + 1/2), fy ^ Ti/{N + AI + 1/2), and 
= Ti/{2N + 1). Specifically, in the limit of large 
squeezing it is predicted that the transverse decay time 
Tx is increased beyond the value of 2Ti owing to the re- 
duced fluctuations in the Q quadrature of the vacuum. 

A simplified schematic of our experiment is shown in 
Figure lb. We realized an effective two level system 
using the ground state and lower energy level of a po- 
lariton formed by a superconducting transmon^^ circuit 
resonantly coupled to the TEioi mode of a 3D supercon- 
ducting cavity.^^ The transition frequency of the effec- 
tive qubit was ujJ2t: = 5.8989 GHz with Ti = 0.65(2) 
/iS set by deliberate coupling to the 50 f2 environment. 
In the supplementary information we show in detail that 
the radiative interaction of the polariton with squeezed 
vacuum is that of an idealized atom interacting directly 
with a squeezed reservoir. Squeezed vacuum was gener- 
ated by pumping a lumped-element Josephson Paramet- 
ric amplifier (LJPA) with two tones at frequencies wi 
and UJ2 that were evenly spaced about the qubit transi- 
tion frequency'^'^ , and satisfied wq = (wi -I- ui2)/2 = Wq. 
The bandwidth of the squeezing was 13 MHz, sufficient 
to fulfill the large bandwidth assumption based on the 
radiative linewidth of the qubit 7/277 = 240 kHz. The 
output of the amplifier was connected with coaxial ca- 
bles to the strongly coupled port of the superconducting 
cavity. 




FIG. 2: Transverse decay into squeezed vacuum, a The 
Ramsey measurement as a function of angle consisted of a 
first 7r/2 rotation about the — x cos 4>-\-y sin (j) axis, to prepare 
the in the state |7r/2,0), followed by a second 7r/2 rotation 
about the —xsinui^adt — ycosuj^nodt axis applied at variable 
time t. The 2D plot displays (az) as a function of t and 
(j) which is characterized by sinusoidal decay with a uniform 
decay constant T2* and phase (j>. b The transverse decay into 
squeezed vacuum was measured by turning the pump for the 
LJPA on interim the qubit pulses. The 2D plot indicates that 
after rapid decay of coherence along the zby axes, the resulting 
coherence along the ±4 axes decays with time constant > 
T2 ■ c The Ramsey measurement for the qubit prepared along 
the —y {(j) = tt) axis with the squeezing off. d,e The Ramsey 
measurement in the presence of squeezed vacuum for the qubit 
prepared along the —y {4> = tt), and +x {(f) = 7r/2) axes. 



To demonstrate the effect of squeezed vacuum on the 
transverse decay of the qubit we conducted Ramsey mea- 
surements at different angles along the equator of the 
Bloch sphere. The Ramsey measurements consisted of an 
initial 7r/2 rotation about the — x cos + ?; sin (/> axis, fol- 
lowed by a second 7r/2 rotation about the — a; sinwmod^ — 
y coswinod^ axis applied at variable time t. Modulation of 
the rotation angle of the second 7r/2 pulse at frequency 
Wijiod results in oscillatory Ramsey fringes without de- 
tuning. Figure 2a displays {uz) as a function of time 
and angle with the squeezing turned off; (ctz) exhibits 
exponentially damped, sinusoidal oscillations at angular 
frequency Wmod and phase (/), with a uniform decay time 



TI = 1.08(4) ^s. 



The reduction of the T2 time from 



2Ti indicates the presence of a small amount of pure de- 
phasing characterized by a time scale = 6.6(5) /zs. 
Figure 2b displays the results of the Ramsey measure- 
ment when the LJPA pump was turned on to generate 
squeezed vacuum for the variable duration between the 
first and second 7r/2 pulses. The power gain of the am- 
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plificr was 4 dB. The transverse decay in the presence of 
squeezed vacmim reveals two timescales, = 1.67 /zs, 
and Ty = 0.28 /is. that describe the exponential decay of 
coherence when the qubit was prepared along the ±x and 
±y axes respectively. Subtracting the pure dephasing 
from the measured timescales gives the radiative trans- 
verse decay times, = 2.2 /xs, and Ty = 0.29 /its. The 
interaction with squeezed vacuum both enhances decay 
along the y axis due to the increased fluctuations in the 
/ quadrature of the field, and suppresses decay along the 
X axis due to the reduced fluctuations in Q. 

The radiative decay dynamics in the presence of 
squeezed vacuum can be presented as a trajectory of the 
Bloch vector. To illustrate this point, we prepared the 
qubit in an initial state \i) with a O.GTtt rotation about 
the -xsincp + ycoscp axis with (p = 0.327r. After this 
preparation, the pump of the LJPA was turned on for 
a variable period of time to generate squeezed vacuum. 
After this variable duration, we tomographically recon- 
structed the qubit state using 7r/2 rotations around the 
X and y axes followed by state readout in the CTz basis to 
determine the Bloch vector components (uy) and {(Jx), 
or no rotation to determine (ct^) . The trajectory of the 
Bloch vector, displayed in Figure 3, follows the expected 
decay dynamics based on (1) with fast decay along the 
y and z axes and by slow relaxation along x. The final 
state of the qubit is described by {(Tx) 0, {(Ty) — > 0.07, 
(cz) — )■ 0.36. The steady state value of (ctz) is consistent 
with a bath of = 0.88 photons that characterize aver- 
age photon occupation of the squeezed state. The rem- 
nant coherence along the y axis is the result of a small 
coherent component of the squeezed state. This coherent 
drive, characterized by a Rabi frequency JIr <C in 
combination with the radiative decay of the qubit results 
in a steady state coherence, {'^y)ss c< ^rTi. Based on 
our measurements, IIr ~ 27r x 10 kHz, consistent with a 
65 dB on/off ratio of our qubit manipulation pulses. 

Because the qubit's decay dynamics are sensitive to 
altered vacuum fluctuations, they can be used to probe 
squeezed states of light. Previously, noise and cor- 
relation measurements have been used to characterize 
the squeezed states generated by microwave parametric 
amplifiers. ^^■^''"^^ Similarly, qubits have been used to to- 
mographically reconstruct localized nonclassical states of 
light. ^^'^^ Here we use the qubit's decay dynamics to to- 
mographically reconstruct, to second order, the Wigner 
distribution for the itinerant squeezed state generated by 
the LJPA. From T^, the measured decay constant of (cr^) 
and we determine N = 0.88 and M = 1.08, from 
which the Gaussian variances aj = (7V + M+ 1/2)/2 and 
(Tq = (A^ — M + l/2)/2 are calculated. Figure 3e displays 
the reconstructed Wigner distribution for the squeezed 
mode at frequency uio- 

In Figure 4a, we display the effective decay constants 
for different values of the detuning 6 = ojq — uiq between 
the center frequency of the pump and the qubit. and 
Ty depend strongly on the detuning near resonance high- 
lighting the "smoking gun" evidence of interaction with 



squeezed vacuum where T^ > 2Ti . When the detuning is 
large, the squeezing axis rotates rapidly with respect to 
the qubit axis and the decay times approach a constant 
value of 2Ti/{2N + 1). The solid black and red lines 
indicate the expected dependence of T^ and Ty on S as 
discussed in the supplementary information. 

In Figure 4b we display the decay constants measured 
for different bias conditions of the LJPA obtained by 
changing the power of the pump tones. The transverse 
decay rates were measured as depicted in Figure 2. As ex- 
pected, larger gain of the LJPA results in larger amounts 
of sqiieezing with an associated increase of T^ and de- 
crease of Ty and T^. Figure 4c displays M — N versus 
N. The reduction of M from its maximum allowed value, 
shown as a dashed line, may be attributed to two possi- 
ble sources; losses in the microwave components between 
the LJPA and the qubit, and non-ideal performance of 
the LJPA. If we assume that the LJPA produces an ideal 
squeezed state, with M = y'A(A -I- 1), then the degra- 
dation can be accounted for by an attenuation of the 
squeezed vacuum from the LJPA by a factor of 77 = 0.5. 
Attenuation degrades the squeezed vacuum by absorbing 
correlated photons thereby making the quadrature fluc- 
tuations tend toward the normal vacuum fluctuations. 
This level of attenuation is consistent with the antic- 
ipated insertion loss between the LJPA and the qubit 
due to the microwave components we used. At values of 
N > 1, however, it appears that the performance of the 
LJPA may become non-ideal as indicated by the slight 
reduction of M - iV for A?" > 1. 

Our results demonstrate the ability to alter the vac- 
uum environment of a two level atom to a degree that 
has so far been elusive in atomic and molecular systems, 
allowing the direct study of a long sought physical phe- 
nomenon of the light-matter interaction. Our system 
also demonstrates the strength of using superconduct- 
ing artificial atoms as sensitive detectors of the quantum 
states of the electromagnetic field. Future studies with 
squeezed light and superconducting qubits may enhance 
the fidelity of quantum gates, enable the generation of 
multi-qubit entanglement^" and enable the study of non- 
Markovian quantum baths. 

Methods Summary 

Due to the finite temperature of the 50 fl environment and 
other sources of noise, a small average number of photons, 
iVth, are expected to contaminate the vacuum environment 
of the qubit. This bath of thermal photons both reduces 
the measured energy decay time Ti from its intrinsic value 
by a factor of l/(2A'th -I- 1), and increases the equilibrium 
excited state population. We determined the equilibrium ex- 
cited state population to be 1.8% using a Rabi population 
measurement^^, allowing us to place a pessimistic limit on 
the number of thermal photons that characterize our vacuum 
environment of Nth < 0.019 and thus the intrinsic radiative 
decay time Ti < 0.67 fis. While a small effect, these thermal 
photons were included in our determination of N and M. 

The polariton qubit was composed of a transmon circuit 
with charging energy Ec/h = 208 MHz and Josephson en- 
orgy Ej/h = 23.27 GHz coupled to a 3D aluminum cav- 
ity with resonance frequency L0c/2n = 6.0456 GHz at rate 
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FIG. 3: Radiative decay dynamics in squeezed vacuum, a-d Quantum state tomography shows the evolution of the Bloch vector 
which was initiaUzed at \i} — |0.677r, O.SSvr) with tomographic measurements equally spaced between and 3 ^is. The dynamics 
are characterized by fast decay along y and z with slow decay along x. e From the radiative decay rates we tomographically 
reconstruct the Wigner quasiprobability distribution of the itinerant squeezed vacuum mode at oJo- The inset shows the 
Gaussian half width of the squeezed and vacuum states as solid and dashed lines respectively. 




FIG. 4: Dependence of the transverse and longitudinal decay times on LJPA detuning and bias, a Effective decay constants vs 
the detuning of the center frequency of the LJPA from the qubit, 5 = luq — iOq. and Ty show a dependence on the detuning 
of the squeezing from the qubit transition frequency, reaching their maximum and minimum values on resonance. The vacuum 
decay limit, 2Ti is shown for comparison. The solid black and red lines are the theoretical dependence of and Ty on the 
detuning. The upper cartoon indicates how the detuning of the squeezing causes the squeezing ellipse to rotate relative to the 
qubit coordinates, b Measured values for Tx, Ty, T^ for increasing gain of the LJPA vs A''. Error bars indicate the standard 
error of the mean based on 10 successive measurements. The upper inset indicates how the aspect ratio of the squeezed state 
changes for increasing A^. c M — A is plotted versus N . The dashed line indicates a minimum uncertainty squeezed state that 
is expected for ideal squeezing. The solid line indicates the expected dependence for a quantum efficiency of 77 = 0.5. The 
gray region indicates values of M and A that are forbidden and the red region indicate values of M and N that correspond to 
classical states of light. 



g/2-K = 126 MHz. The cavity was equipped with two ports; 
one strongly coupled port that limited the quality factor to 
Q = 1.1 X 10* and another weakly coupled port. The qubit 
was enclosed in successive layers of superconducting and mag- 
netic shields and anchored to the mixing chamber stage of a 
dilution refrigerator with a base temperature of 20 mK. State 
readout was performed using the Jaynes-Cummings nonlin- 
earity technique'^'^ by driving the weakly coupled port of the 
cavity with a strong tone at 6.0467 GHz and integrating the 



first 200 ns of transmitted signal. 

The LJPA was composed of a two junction SQUID formed 
of 1 ^A Josephson junctions shunted with 1 pF of capaci- 
tance and isolated from the input ports of a 180° hybrid with 
interdigitated capacitors that resulted in a quality factor of 
Qljpa = 100. The LJPA was flux biased to have a low power 
resonance at 5.897 GHz. The differential port of the hybrid 
was connected to the strongly coupled port of the qubit with 
coaxial lines via two circulators and a -20 dB coupler that 
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allowed the injection of the qubit manipulation pulses. 

A single microwave source was used to generate the qubit 
preparation, tomography, and LJPA pump pulses. The LJPA 
pump was obtained by driving an IQ mixer with a tone at 
540 MHz, and adjusting the dc offsets to null the carrier. 
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I. THE MASTER EQUATION FOR THE TRANSMON-CAVITY SYSTEM 

In the main text, we state that the irreversible dynamics of the effective two-level system (TLS) defined by the 
ground state and lower polariton of the transmon-cavity system are the same as for an atom interacting with a 
broadband squeezed vacuum. In this section, we give the theoretical reasoning behind this statement. We derive the 
effective master equation for the transmon-cavity system and compare it to previous predictions for an atomic system 
interacting directly with the squeezed vacuum. This problem was treated before in the "bad cavity" limit where the 
cavity decay was taken to be the fastest process in the system [1]. Here, this assumption does not hold and we follow 
a different derivation. 

The system is a transmon circuit strongly coupled to a superconducting cavity. The cavity is also coupled to the 
line modes in which the squeezed reservoir is realized. The chosen system parameters allow for several theoretical sim- 
plifications. The squeezing bandwidth [BW = 13 MHz) was 50 times the radiative qiibit lincwidth (7/27r = 240 kHz). 
From the point of view of the qubit, the squeezing is well approximated as having infinite bandwidth. However, since 
the polariton transitions are well separated with respect to the squeezing bandwidth, with 255 MHz between the 
\g) — )■ |— ) and \g) — >■ |-|-) transitions, the squeezed vacuum only addresses one pair of levels, an approximation that 
we make in Section 2. The several-GHz qubit transition frequencies allow us to use the rotating wave approximation. 
We additionally assume the Born approximation, which terminates perturbation theory at the second order, and the 
Markov approximation, which assumes that the squeezed reservoir has no memory (its correlation time approaches 
zero). 

The Hamiltonian for this system is 

U=U^(^ '^kblbk + ^ hfe (bl + bk) (at + a) (1) 

k k 

where is the transmon-cavity Hamiltonian, ojk is the frequency corresponding to the k^^ traveling wave line 

mode, bj. (bk) is the creation (annihilation) operator for the fc**^ line mode and (a) is the creation (annihilation) 
operator for the cavity photon, h^ are the real matrix elements of the transitions between the k^^ line mode and 
the cavity mode. The transmon-cavity Hamiltonian Ti^'-^ is a generalized Jaynes-Cummings Hamiltonian for the 
transmon levels and the cavity mode [2] . 

To express the system evolution in terms of the line mode correlations, we first diagonalize the evolution to represent 
it in a joint polariton and line mode basis [3]. We apply the transformation U that diagonalizes the closed system 
part of the Hamiltonian H'^'-^ to H 

H' = no + J2hk (bl + bk) U (at + a) (2) 

k 

= no + Y.hk (bl + bk) (^^.- N) (il + li) (^1) • (3) 

k ij>* 

Here, the free evolution is Ho = J2i^I'" (*l + J2k'^kblbk and we express the transformed cavity operators in 
terms of ladder operators between the polariton levels |i) , |j) in the diagonalized system U [a + a^) W = + a = 
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Si j>i i^ij I*) ^-^1 ^ij This expansion relies the herniiticity of + a to allow the restriction j > i and since 

+ a changes the parity of the state, An = [3]. The energy of the i^^ polariton is denoted ef 
Now, we move to the interaction picture with respect to Ho where we can apply perturbation theory. The coupling 

term V = Y]^. hfe (bl + 6fe) {Aij \i) {j\ + A*ij \j) is written as 



(4) 



Here, we define a modified line mode photon annihilation operator Bt = X^jt ^kC- "^''°*&/c and the polarition jump 



operator at = Ylii j^i-^ij^ K) 01- We have used the rotating wave approximation to remove the terms -B^crJ and 
Btut. The energy difference between polarition levels is written as Aj^ — Cj — 

In perturbation theory, the system evolution is = i [V/,Pj]. Applying this twice, we obtain 

V/, /* dTi [Vi, = f dr {Vipl-Vl - Vlvip[ + h.c.) 

^0 J ^0 



(5) 



Taking the partial trace over the line modes, we get the evolution of xc ~ ^rs [p^] . To evaluate this, we expand 
using Equation (4) and use the fact that the trace is cyclic. We further assume the Born approximation [5] 

'bIbI' 



= Pxc ® Pb to evaluate terms like Trs Blarp^Blat = (TtP^ ^jr^atTiB B\p^^ ^Bj 

Pt,TC = { (o-rPx.TcC't - CTtCT^p^ Tc) (b\b\^ + {a].pl.^T^c<^l - ala].p[^T^(}j {BtB^} 

+ (4Pr,TC0-t - (Ttalpi^^c) (bIBt) + (o-rPr.TCO-l " (^UTpi,Tc) {^tBl) + h.C.j . 

Consider one of these terms, Jq dr ajarP^'j^c (^BtBiy Expanding the reservoir correlator (^BtBiy. 

(BtBl) = /^hfce-'"'=*6fc^hfe.e'-'=-4,\ 

\ fe k' I 

= ^hfchfe.e-'(-'=*-'=-)TrB [{bibk'+Skk') pi,B 

k,k' 

= ^(iVfe + l)h2e--'=(*--) 

k 

where 6kk' is the Kronecker delta. Incorporating this into the full term 

(' dT<j\arp[,rc{BtB\) = Tdr ^ A*^.e-'^-^*|j)(i| 

X J2 Ai'i'e'^^v-|i')0-V^_^<,^(iVfc + l)h2e--'=(*--) 



(6) 



(7) 
(8) 
(9) 

(10) 



i,j>i k 



(11) 



where = {h\hk)- We will neglect fast-oscillating terms which have little effect on the decay rates by effectively 
taking j = j'. Other terms will give rise to the squeezing correlation function {hkhk') = Mkiu^^+^oy ,2LJa [6]- The sum in 
this term (11) is an expression that refers to the line modes. The temporal kernel in the integral (10) can be written 
in frequency space as 

/oo 
dujk D (wfe) (iV„, + 1) h2^e--('^'=+^*^-)(*-") (12) 
-°° 

where D (•) is the density of states at the given frequency. Since we are interested in the case of a broadband reservoirs 

Ni^^ has a bandwidth large compared to the typical decay rates and therefore the function K{t—T) decays on a timescalc 
of the order of 1/ BW. In the limit of infinite bandwidth it acts similarly to a Dirac delta function effectively restricting 



3 



the integral to t « This behavior justifies applying the usual set of Born-Markov approximations [7], under which 
the term (10) becomes 

E l^'il' S^jSr.pl^c [N^^^ + 1) (13) 

where S^j = \j) {i\ and S^^ = \i) {j\ with j > i and where 7,j = 2-kD (A^j) /i^(Ajj). 

Computing the other terms similarly and transforming back from the interaction picture, we find 

Pt,TC = Pt = -i[^o,Pt] (14) 

i,j>i 

+ {2Sr-.pSr. - Sr.Sr.p - pSr.Sr^) ^Ml^^.e^H-o+A.,)* 
+ {2Sr.pS± - Spr.p - pS±Sr^) |A,,f (^a., + l) 
+ {^S±pSr. - Sr.S±p - pSr.St^ |A,,f J 

where wq is the center frequency of the squeezed reservoir. This equation is the multi-level generalization of the 
two-level master equation derived by Gardiner, Eq. (10) of Ref. [4] [8]. 



II. DETUNING DEPENDENCE OF POLARITON COHERENCE DECAY 

When we restrict the master equation (14) to only two levels, we obtain a mathematically equivalent master 
equation to Ref. [4]. Physically, the two lowest levels i = 0, 1 refer to the ground state \g) and the dressed state 
|— ), respectively. In the experiment, the squeezed vacuum bandwidth is centered around the \g) — )■ |— ) transition 
frequency. Since all other transitions are well outside this bandwidth, we can limit ourselves to the two-level dissipative 
dynamics and directly apply Gardiner's theory to the ground state and lower polariton of the circuit QED system 
explored experimentally in the main text. 

In the absence of a classical coherent drive, the decay of atomic polarizations obey on average 



&+ = -TN<J+ + TMe^'^'a_ (15) 
a-=-VNa-+Tl,e-^'^'cj+ (16) 

where (T± = {S^i), Vn = l(N + 1/2) + 7^, Tm = iM, 7 = |Aoip7oi, and 6 = ojq + Aqi = + Wq. 7^ is a 
phenomenological pure dephasing that was measured independently. 

The decay dynamics of the atomic coherences ax,y can be found by solving Equations (15) and (16) and then 
transforming to a frame rotating with angular frequency 5, i.e. (T±{t) = e^^^^a±(t). These dynamics will have two 
decay timescales Tn ± vITmF^^ = l{N + 1/2) + 7^ ± s/^f\M^f^^ unless 5 = Q ov\8\> ^\M\. Detuning also 
leads to an effective time-dependent rotation of the squeezing ellipse. This limits the influence of squeezing on decay 
to the range \5\ < j\M\ even in the ideal case of infinite bandwidth squeezing. 

To interpolate between the limits of large detuning {\d\ > 7|M|) and 6 = where the decay is described by 
a single exponential time constant. Figure 4a in the main text presents the effective decay constants that were 
obtained by fitting the Ramsey measurements to a single exponential decay at each detuning. For comparison to 
the theoretical prediction, the same fitting routine was applied to theoretical decay curves given by Equations (15) 
and (16) by transforming into the frame rotating with angular frequency S + Wmod, with u^od/'^'^ = 5 MHz and 
70 = 27r X 2.4 x 10^^ s^^. The decay constants that result from this fit are plotted as the theoretical curve in Figure 
4a. For a direct comparison to theory, we present experimental data and the predicted Ramsey measurement results 
versus detuning in Figure SI. 
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